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Abstract 

We calculated the contribution of the photon hadronic component in the charmed 
structure function of the leptoproduction. The contribution comes from scattering of the 
c quarks of the virtual photon on the quarks and gluons of the proton. Comparison of our 
calculations with the measurements of the charm production in /u + p-scattering by EMC 
shows, that contribution of resolved photon can explain the excess of EMC data over the 
predictions of the photon-gluon fusion model at large momentum transfers. Thus, one 
does not need to use a non-perturbative admixture of the charmed quarks in the proton 
wave function ("intrinsic charm") to describe such an excess in EMC data. 
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1 Introduction 



Investigations of heavy flavours production at high energies remains to provide a very important 
tool for qualitative and quantitative test of the QCD and for study of the internal hadron 
structure. Despite the impressive experimental and theoretical results obtained during last 
decade due to, mainly, experiments at the electron-proton collider HERA, there still exist 
uncertainties in the interpretation of data on the charm production in hadron-hadron and 
lepton-hadron collisions. These uncertainties relate to the details of the quark structure of 
hadrons, in first turn the proton. One of these problems is the problem of "intrinsic charm" 
(IC) of the proton [1], i.e., the question about the presence in the wave function of the hadron 
a visible 1%) non-perturbative admixture of the charmed quarks with hard, "valence- like" 
distribution over their longitudinal momenta. Investigations at HERA have been performed at 
small values of x, where the main contribution comes from the photon-gluon fusion mechanism 
(PGF). But small values of the variable x describe the space scales much greater, than the 
proton size and, therefore, the study of those kinematical regions gives information about the 
structure of the QCD vacuum itself, rather than about the internal structure of the proton. 
To understand the internal structure of the proton it is necessary to study the proton charmed 
structure function at large values of x. But the geometries of ZEUS and HI set-ups do not 
allow to study efficiently the charm production processes in the very forward region at relatively 
small momentum transfers, where the cross-section of the ep scattering is maximal. The planned 
HERA upgrade will increase few times the collider luminosity and, perhaps, this increase will 
allow to collect sufficient statistics at large momentum transfers for the experimental study of 
the charm production in the range of large values of x. But at present this problem is not well 
investigated. 

The theoretical consideration of the "intrinsic charm" problem is based on the following 
simple points. The wave function of the proton can be expanded over the colourless eigenstates 
of free Hamiltonian \uud >, \uudg >, \uudqq >, .... During sufficiently long time the proton 
can consist of Fock states of an arbitrary complexity, including a pair of charmed quarks. In the 
proton rest frame the life-time of such a fluctuation r has the order of the nuclear time ~ Rh, 
where Rh is the hadron size. Charmed quarks are heavy objects and their life-time is much less 
than that of light partons. Therefore, in average, the admixture of the heavy quark pairs in the 
proton wave function is to be small, ~ (m,/mg) 2 , that is, by factor ~ 10 2 less than life-time 
of a fluctuation containing light partons only. Additionally, because of c quarks are massive 
it turns out that their momentum distributions are to be much harder than the distribution 
of the light sea partons [1, 2]. The existence of the fluctuation, containing heavy quarks is a 
natural consequence of the field theory. But the important guess is that quantum fluctuations 
of the hadron wave function containing charmed quarks appear due to self-interaction of the 
colour field, ensuring the non-perturbative origin of the intrinsic charm distribution. In this 
case the structure of the hadron Fock states with charm can be considered independently of 
hard interactions. 

Estimating the level of the charmed quarks admixture in the proton, they usually suppose 
that experimental indications on the intrinsic charm exist in [ip collisions [3] and in hadronic 
collisions data [4, 5] , where deviations from the photon-gluon and parton-parton fusion models 
predictions have been observed. But the experimental situation is sufficiently uncertain. The 
original work [1] was motivated by ISR data [4] , where a very large yield of the charm had been 
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obtained, at least, order of magnitude larger than it was expected in the model of the parton- 
parton fusion. A critical comparison of the ISR results against the contemporary fixed-target 
data can be found in the review [6]. Later experiments (though carried out at lower energies), 
did not contradict so dramatically to the theoretical expectations (see, e.g., [7]). In paper [8] an 
attempt has been taken to co-ordinate the ISR data with the parton-parton fusion predictions. 
There has been considered the process of the "charm excitation", based on the guess about 0.5% 
admixture of the charmed quarks in the proton, which admixture leads to the hard scattering 
of the charmed quarks. Due to large cross section of the hard scattering (strictly speaking, this 
cross section diverges) and to its stronger dependence on reaction energy, it was shown [8] , that 
ISR data can be described in framework of the charm excitation model. At the same time this 
model does not predict visible contribution in the charm yield at lower energies. 

Experimental data on longitudinal distributions of charmed particles and on asymmetry of 
their production can be reproduced with acceptable accuracy without the "intrinsic charm" 
hypothesis, only by reasonable variation of parameters of string hadronization model [9, 10]. 

The interpretation of the data on J/ip particle production is also contradictory. In the 
experiment NA3 [5] in pion-nuclear collisions at energy 280 GeV definite indications on the 
presence of an additional mechanism of hard J / ip production were obtained. This additional 
contribution at the level of 20% of the total of J ftp production cross section is well described 
by the modified model of the intrinsic charm [2], but the guess about so large contribution does 
not look very realistic. It was noted in papers [11, 12], that intrinsic charm can also explain a 
larger than it was expected yield of fast correlated pairs of J / ip particles, as well as anomalous 
polarization of J/V's, probably observed in NA3 experiment. At the same time predictions 
of the intrinsic charm model directly contradict to data on J/ip production in proton- nuclear 
collisions at energy 800 GeV [13] (see [14]). 

Note, that the admixture of the charmed quarks in the proton at level of 1% leads to total 
cross section of charm production in nucleon collisions of order of few hundreds microbarn in 
conflict to the bulk of the experimental data on open charm production in hadronic interactions. 

Doing any theoretical analysis of the hadronic data it is necessary to keep in mind, that 
the main contribution in the charm production in hadronic collisions is led by the gluon-gluon 
fusion gg — > cc. At the same time the direct measurements of the gluon distributions are 
impossible. To extract the gluon distributions from data they use either charm production 
calculations based on the gluon-gluon and photon-gluon mechanisms, which contributes in the 
small xf range, or sufficiently involved methods based on the guess that the main contribution 
in the QCD evolution is driven by the gluons. Due to this fact the extraction of an additional 
mechanism of charm production at level of (0.5 — 1)% can not be reliable. 

In fact, the only experimental data, which allow to perform sufficiently reliable check of 
the presence of the intrinsic charm in the proton are the EMC measurements of the charm 
photoproduction in fx + p collisions at energy 250 GeV. In early papers on the EMC data analysis 
they used old, harder parametrizations of the gluon distributions in the proton. In particular, 
in the first analysis carried out by EMC the scaling gluon distribution G(x) ~ (1 — x) 5 has been 
used, which gave the deviation from the photon-gluon mechanism predictions at level of about 
0.3% with sufficiently large error of this value. Modern and, as a rule, softer parametrizations 
of the gluon distribution, lead, naturally, to larger contribution of the intrinsic charm in the 
EMC data [15, 16, 2]. Thus, the intrinsic charm contribution, extracted from analysis of the 
EMC data varies from 0.3% [3, 17] to, approximately, 0.9% [16, 2]. In paper [16], where the 
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most careful calculations of the contribution of the photon-gluon fusion have been done, it 
was shown that at small energy transfers the intrinsic charm contribution is absent within the 
experimental errors. 

Due to existence of the hadronic structure of photon [18, 19, 20] the total charm production 
cross section in the leptoproduction should also include the interactions of the partons from the 
virtual photon with the proton. Such a contribution of the hard scattering of c quarks from the 
virtual J /ip meson calculated in framework of the vector dominance model in [21], allowed to 
reproduce with a good accuracy the results of ZEUS experiment on D* mesons production at 
large transverse momenta. Despite the fact that the scattering processes are to the second order 
by a s , they can produce quite visible effects due to divergence of the parton-parton scattering 
cross sections at small angles. Therefore, it is natural to expect, that the hadronic component 
of the resolved photon can also give a visible contribution in the charmed structure function, 
measured by EMC in the range Q 2 > 40 GeV 2 , where the indications on the deviation of the 
experimental data from the photon-gluon predictions exist. 

In present work we have calculated the contribution of the hadronic component of the 
photon and the contribution of the standard mechanism of the photon-gluon fusion in the 
charmed structure function, measured in EMC experiment. The calculations performed in 
present work show, that it is not necessary to consider the intrinsic charm admixture in the 
proton to interpret the EMC data. 



2 The basic outline of the model 



2.1 Charm in the hadronic component of the photon 

In the deeply inelastic processes with hadrons the measured physical quantity is the hadron 
structure function F 2 (x, Q 2 ). The charm production cross section in electroproduction hp 



expressed through the structure function F 2 c (:r, Q 2 ) is: 

27ra 2 



da(lp — ► c) = 



xQ l 



[l + (l-y) 2 ] F^(x,Q 2 )dxdQ 2 , 



c, 



(1) 



and in the equivalent photon approximation [22] it can be expressed through the cross section 
of the 7p interaction a(jp — > c) and the flux of the equivalent photons n y (x, Q 2 ) as: 



da(lp — ► c) = dn 1 (x,Q 2 ) a(jp —> c) 



dn 7 (x,Q 2 ) = § 1 + (1-y) 



dxdQ 2 
x Q 2 



(2) 



That is, the charmed structure function is connected to the charm photoproduction cross 
section as the following: 



F 2 c (x,Q 2 ) = 



4tt 2 



a 



<j(7P -> c) 



(3) 
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The hadronic structure of the photon can also be described by the photon structure function 
F^x^Q 2 ), depending on Bjorken variable x and on the momentum transfer squared Q 2 [18]. 
The structure function F^(x, Q 2 ) is usually subdivided in two terms - perturbative (so called 
"anomalous") and non-perturbative ("hadronic") parts. Such a subdivision of F^(x, Q 2 ) is valid 
not only within the framework of the naive quark parton model (QPM), but also in the next-to- 
leading order of QCD [18, 19]. The perturbative part is led by direct interaction 7 — > qq and 
can be completely calculated in the framework of the perturbation theory. The contribution of 
this part in QPM is proportional to In Q 2 and dominates at Q 2 — > 00. The non-perturbative 
part of F^x^Q 2 ) is completely analogous to the ordinary structure function of hadrons. In 
the framework of the naive QPM this part does not depend on Q 2 . The QCD evolution of 
F2 differs from the evolution of the hadron structure function only due to the existence of the 
perturbative part in the F 2 7 thus being described by the inhomogenious DGLAP equations. 
To calculate the QCD evolution of the non-perturbative part, similarly to the hadron case, it 
is necessary to set the initial conditions at some value of momentum transfer Qq. As a rule, 
choosing the initial conditions for the parton distributions in the photon at Qq ~ 1 GeV 2 
they suppose, that the hadronic part of F]{x, Qq) can be derived from the vector dominance 
model (VDM), postulating the transition 7 — > p° to describe the interaction of the photon with 
hadrons. The parton distribution in p° meson are obtained using the parton distributions in 
pions known from the experiment. For better description of the experimental data the simple 
VDM has been extended by the obvious way by inclusion of iso-scalar vector mesons u, <fi. The 
extention of the VMD by inclusion of heavier mesons from J/ip family is not, normally, used 
because of large masses of J/ip mesons and the charm content in the photon is described by 
the evolution of light partons and by the perturbative term 7 — > cc [20, 23]. But, as it has 
been shown in [20] , in the range Q 2 < 50 GeV 2 charm must not be considered as a light quark 
and the evolution equations must be solved for three light flavours only. This is because of the 
fact that in the threshold region Q 2 ~ Am 2 the charm production must be described by the 
Bethe- Heftier process 7 7 — > cc and by its analogue, the photon-gluon fusion process 7 g — > cc 
in case of leptoproduction. On the other hand, in paper [24] from comparison of the QCD 
evolution predictions for F^ to experimental data the conclusion has been made that even for 
evolution of light quarks the perturbative region in case of the photon starts at Q\ ~ 5 GeV 2 , 
not at Qq ~ 1 GeV 2 . Thus, in present work we use the generalized model of vector dominance 
(GVMD), taking into account also non-perturbative transitions of the photon to J ftp mesons, 
which become possible beyond the J/tp production threshold. 

Then, it is suitable to represent the wave function of the hadronic component of the photon 
as following: 



V V JV q 

(c) 

where terms |7 per t > and \T pf ! rt > take into account the perturbative contributions of the 
processes 7 — > qq and 7 — > cc, respectively, and vector mesons V = p, u, 0, J/ip describe 
the non-perturbative contribution in the hadronic component of the photon. The constants, 
characterizing the contributions of vector mesons are [23]: 
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/ p 2 /4tt « 2.20, /2 /47r ~ 23.6, 

(5) 

/2/47T « 18.4, /|/47T « 11.5. 

At presence of the hadronic structure of photon, the interactions of its hadronic component 
can be described in analogue to the ordinary hadrons as interactions of the partons from the 
photon with the colliding particle. In the framework of the parton model the cross section of any 
process is a sum of the cross sections of all elementary subprocesses, weighted by distribution 
functions of partons, participating in these subprocesses. If charmed quarks exist in the proton 
(" intrinsic charm" ) the leading over strong interaction constant process at moderate momentum 
transfers is the scattering of the lepton on charmed quark of the proton Ic — > Ic. The photon- 
gluon fusion process 7 g — > cc has the next order in a s , but its contribution dominates due to 
the fact that gluons carry almost half of the proton momentum. Thus, one can expect that 
the contribution in the charm photoproduction comes from four processes - the absorption of 
the virtual point-like photon by c quark of the proton, the photon-gluon fusion, the scattering 
of the c quark from virtual photon on light partons of the proton and the scattering of light 
partons from the resolved photon on c quarks of the proton. Therefore, the total cross section 
of the charm production can be written as: 

a( w ) = a IC + a PGF + 0* + a® + a® (6) 

Here, a IC defines the absorption of the virtual photon by c quark of the proton; the cross 
section a PGF describes the photon-gluon process 75 — > cc. The second pair of terms in (6) 
comes from scattering of c quarks from the photon on the partons of the proton and a® is led 
by scattering of light partons from the photon on c quark of the proton. The light partons of the 
photon can originate from the non-perturbative part (light vector mesons p, u, 0), as well as 
from perturbative part Yi ver t >• It follows from (5), that the contribution of J/ip meson in the 
hadronic component of the photon is about 16% of the contribution of the light vector mesons. 
At the same time the contribution of scattering of light partons from photon on charmed 
quarks from the proton a® is proportional to the level of admixture of the charmed quarks in 
the proton, N IC < 0.01 << 1, and can be neglected. The term describes the contribution 
of the scattering of the "perturbative" c quark from the photon on partons of the proton. We 
denoted as "perturbative" those c quarks of the photon, which appear either due to anomalous 
part of the photon hadronic component or to QCD evolution of original light partons. Thus, 
the distribution of such c quarks can be obtained in the framework of the perturbation theory 
[20, 25] and, as quoted above, this contribution is limited by the range Q 2 > 50 GeV 2 . Using 
the parametrization [26] for the distribution of the "perturbative" c quarks 1 , we have verified, 
that the contribution of the perturbative term can be safely neglected in the kinematical 
range considered in present work, Q 2 < 80 GeV 2 . 

The contributions in the charm production come also from the processes of parton-parton 
fusion gg — > cc and qq — > cc, completely analogous to the processes, describing the charm 
production in hadronic collisions. But, these processes have finite total cross sections and their 
contribution in the charm yield is small in comparison to the photon-gluon fusion process. Thus, 

technically, we used the respective subroutines from the package PYTHIA 6.1 [27]. 
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we also drop them. Quite different situation appears for scattering processes of c quarks from 
J/ip meson on light partons of the proton. In the pQCD the cross sections of these processes 
are divergent as Therefore their contributions can be significant. As the result of the 

above consideration we keep in (6) only terms a IC , a PGF and (7^. Their contributions in the 
charmed structure function will be calculated further. The diagrams of the processes, taken 
into account in present work are depicted in Fig. 1. 



2.2 The c quark distribution in virtual J/ip meson 

To find the non-perturbative distributions of c- and c quarks in the virtual J / ip meson we can 
follow the approach developed in [2] . This approach is based on the statistical model [28] and 
on the work [1], allowing to take into account the masses of heavy quarks in the framework of 
the non-covariant perturbation theory (NPT). This approach allows to find the probability of 
creation of the Fock state with heavy quarks. The basic point is the fact, that we can use NPT to 
derive the necessary distribution. Feynman diagram is a sum over all time-ordered diagrams of 
the NPT. The contributions of time-disordered diagrams are proportional to the inverse powers 
of large hadron momentum P^. Thus, any consideration in the framework of the parton model 
is being carried out in the "infinite momentum reference frame", i.e., in the reference frame, 
where the hadron momentum is much larger than any of characteristic mass parameters. At 
high energies of the colliding particles the momentum of the hadron is sufficiently large, which 
ensures the applicability of NPT. Only in this case the virtual configuration of Fock states, in 
which the proton fluctuates, can be "frozen" for the time of the interaction. If we consider the 
processes led by the transition of the photon into a hadron (vector meson), the life-time of the 
virtual hadronic fluctuation (given Fock state) is At ~ 1/AE ~ 2P h / (M 2 — m 2 ) (P h - the 
hadron momentum, m - hadron mass and M - the mass of the fluctuation) and at high energy 
can be large enough even for large values of the fluctuation mass M. 

Expression (6) means, in fact, that different terms in (4) do not interfere. This, in turn, 
allows the probabilistic interpretation of the expansion (4). Namely, it is legitimate to assume 
that during the interaction time the photon is in one of the states, entering the expansion (4). 
In this case partons of the photon carry all the momentum of the photon (more correctly, all 
the momentum of the hadronic fluctuation). Due to sharp cut-off over transverse momenta of 
partons, it is sufficient to consider the longitudinal phase space only. We consider a hadron (or 
hadronic fluctuation of the photon in our case) as a statistical system, consisting of m partons, 
carrying in sum the quantum numbers of the hadron and its momentum. In spirit of the parton 
model we suppose an independent creation of each parton with relative momentum £ according 
to the probability density Then, in the framework of NPT the probability to find n-parton 
final state can be written as [2]: 

dW^ , 5^ ( P fin - P h ) 17 ^ , (7) 

where, n is the total number of partons in the considered fluctuation, Ph and Eh are momentum 
and energy of the hadron; Pf in and Ef in are total momentum and energy of the fluctuation, 
respectively, (^-function enforces the conservation of the total 3-momentum. Function pi(£i) 
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describes the probability to create iih parton without conservation of the total momentum 
of the system and its quantum numbers. As an argument of the distribution function in the 
virtual J ftp meson we use the light-cone variable & = = (e + p L )/(E h + P h ), where e, p L 
- energy and longitudinal momentum of a parton. 

Because of the life-time of the fluctuation with a pair of charmed quarks being much less 
than the one for the fluctuation, consisting of light partons only the full quark-gluon cascade 
cannot develop. This leads to the situation, where the newly created fluctuation with charmed 
quarks does not practically contain light partons and the pair of heavy quarks carries all the 
momentum of the virtual J/if). Therefore, the probability to observe the Fock state with cc 
pair can be written as follows [2]: 



dwa = eAlp c (Q P ft) s(l _ & _ &) ^ d§ ( 

In J/ip meson charmed quarks are the valence quarks. Thus, coming from Regge phe- 
nomenology, one can expect, that at small values of the relative momentum the probability 
density p(£) can be parametrized as: 



Pc(o = p- c (o = r, (9) 

where, a ~ 0.5, as in ordinary hadrons. So, the inclusive distribution of the charmed quark in 
the photon has the form: 



c 7 (o = e a - +i, (io) 

with (3 = a + 1; B(u,v) is beta-function and the charmed quarks distributions are normalized 
on unity, J* d£c 7 (£) = 1. Choosing a = 0.5, one obtains (3 = 1.5. We shall use this value 
in our calculations. The value of f3 — 1.5 is smaller than, e.g., chosen in [21] for description 
of spectra of charmed mesons, measured in experiment ZEUS. But in our case the correct 
form of the distribution of c quarks does not play any significant role since the contribution in 
the structure function is defined mainly by the value of the total cross section of the parton 
scattering. Then, the momentum transfers in EMC experiment Q 2 < 80 GeV 2 practically, are 
completely in the range, where the perturbative consideration of the production and evolution 
of c quarks is invalid [20]. Thus, we do not consider QCD evolution of c quarks in the virtual 
J/ip meson and we shall use scaling expression (10). 



2.3 Charmed structure function in leptoproduction 

Due to the fact that the structure function is up to a common factor simply sum of cross sections 
of the subprocesses, it can be also represented as a sum of contributions of three processes (see 
(6)): 

F 2 c (x,Q 2 ) = F/ c (x,Q 2 ) + F 2 PGF (x,Q 2 ) + F}(x,Q 2 ). (11) 
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In framework of parton model the structure function of the "intrinsic charm" Fl° is con- 
nected to the distribution of charmed quarks in the proton c p (x, Q 2 ) and is equal to: 

Fi c (x,Q 2 ) = 2N Ic e 2 c xc p (x,Q 2 ), (12) 

where, xc p (x,Q 2 ) is normalized on unity, e c = | is the electric charge of c quark and Nic 
defines the level of the admixture of charmed quarks in the proton. 

The structure function -F 2 PGF defines the photon-gluon fusion and is [29]: 

Ff GF ( x ,(f) = f d A G (tQ 2 )f4,Q 2 ), (13) 

where, 



/ 2 (^Q 2 ) = ^etir zU [-§ + 2z(l - z)(2 - \)] 

1 . (14) 

+ [1 - 2z(l - z) + 4Az(l - 3z) - 8AV] ln}±^ \. 

In the above expressions s = Q 2 (1 — z)/z, V c (s) = \Jl — Am 2 /s is the velocity of c quark 
in the (7*7) center of mass system, A = m 2 /Q 2 . 

The contribution of the virtual J/ip meson, Frf , in the charmed structure function is being 
led by the expression (3) and the cross section of the process J/iftp — > c in framework of parton 
model is: 

°* = 2 ^7?E / ^i^cC^te)^)^). (15) 

h i J 

Factor "2" takes into account the scattering both on c and c quarks. Index % runs over all 
partons in the proton, participating in the scattering off the c quark from the photon. Functions 
G v {^i) and c 7 (£ c ) describe the momentum distributions of the ith light parton in the proton 
and c quark in the photon, respectively. The cross section of the scattering of c quark on the 
parton of the proton di(s) depends on the CM energy squared s. It is known that at presence 
of non-perturbative effects, like quark masses, parton model is not a covariant one. Thus when 
performing the integration in (15) it is necessary to choose some reference frame, satisfying the 
conditions of the applicabilty of parton model. In our case the most suitable reference frame 
is the center of mass system of the virtual photon and the proton. If we use the light-cone 
variable £ + the energy squared s can be expressed through invariant variables x and Q 2 and is 
equal to s = m 2 + £ c £j Q 2 /x. 

Formulae for the cross sections of the scattering subprocesses eg — > eg and cq — > cq of 
massive c quarks in the lowest order over a s can be found in [8, 30]. For the completeness 
we give these formulae in the Appendix. Pure kinematically the minimal momentum transfer 
is imin = 0. Because the cross section of the scattering process is divergent at lower limit 
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of the integration over t, the first question is, obviously, about the cut-off value t m i n . At low 
momentum transfers pQCD becomes invalid and the value of the cut-off momentum also cannot 
be calculated theoretically. For the processes with c quarks various choices of \t min \ are being 
used, varying from m 2 /A to m 2 cT , where m 2 T is the transverse mass of c quark squared. We 
shall consider the value of t min as a free parameter, which must be found from experimental 
data. We use the same value of \t m i n \ also in the argument of the strong coupling constant 
OL s {\tmin\) ■> because the main contribution in the total cross section is defined mainly by the 
range t ~ i min . 

Since we consider the process 7p — * cc + X with large momentum transfer, a coalescence of 
one of the final c quarks with the proton remnant is not probable. Thus, according to quantum 
numbers conservation, there must be as minimum three particles in the final state - the proton 
and two D mesons. So, general kinematical restriction for the scattering reaction takes the 
form {mu is mass of D meson, m p is the proton mass). 

^ ^ Q 2 + m\ > (2m D + m v f . (16) 

At presence of the non-zero cut-off t m in there are additional kinematic constrains on thresh- 
old energies V^i/i [85 30] for the subprocess. We also give them in Appendix. 

It is easy to obtain necessary integration limits in (15), considering the kinematics of the 
parton scattering. Taking into account non-zero value of t m i n the integration limits over ^ and 



X 



1 , (Vsth+m c f 

1 + — w — 



x (sth-mp 



2 



< 6 < 1, 

(17) 

< & < I- 



3 Results and conclusion 

We fitted the calculated F£(x, Q 2 ) to EMC data. The charm photoproduction cross section 
<r(7p), measured in [3] has been recalculated to the structure function as described in [2]. 
Comparing the calculations with data we took into account the next-to-leading order QCD 
corrections to the "intrinsic charm" structure function F 2 /c , as well as non-zero masses of the 
proton and charmed quark (see [2, 17, 31]). When calculating the contribution of the photon- 
gluon fusion we used the parametrization MRS(G) for the gluon distribution in the proton from 
PDFLIB package [32]. The parametrization MRS(G) has been obtained with a s in the second 
order and Aqcd = 0.174 GeV. We included mass of the c quark in free parameters to check 
the self-consistency of our calculations. 

Fit of EMC data gave the following values of free parameters: 
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N IC = (0.2 ± 0.2)%, 

(18) 

m c = (1.51 ± 0.03) GeV, 

(19) 

Lin = (-3.0 ± 0.3) GeV 2 , 

at good value of x 2 /NDF = 0.74. The value of cut-off over the momentum transfer squared 
is approximately equal to the transverse c quark mass squared, as it could be expected from 
general considerations. The resulting value m c = 1.51 ± 0.03 GeV agrees with m c = 1.50 GeV, 
used in the parametrization MRS(G), which testifies the self-consistency of our approach. The 
comparison of the model with EMC data [3] is given in Fig.2. 

Note the following. In paper [8] the calculations have been performed for the charm ex- 
citation mechanism in hadron collisions, i.e., there have been considered the processes of the 
scattering of partons on the charmed quarks of the proton (on "intrinsic charmed quarks"). 
The contribution in the charm leptoproduction formed by hadronic component of the photon is 
completely analogous to the charm excitation in hadron collisions, considered in [8]. In paper [8] 
the level of the intrinsic charm in the proton has been chosen equal to 0.5% and the conclusion 
was made that the charm excitation mechanism ensures the needed yield of charmed particles 
at ISR energies due to steeper rise of the scattering cross section with reaction energy, which 
allows to agree ISR data and data obtained at lower energies in fixed target experiments. In 
paper [8] the cut-off momentum transfer has been chosen equal to i min = m 2 /4 pa 0.6 GeV 2 , 
which is considerably less than obtained in present paper. 

The results of present paper give significantly less admixture of the intrinsic charm in the 
proton and much larger value of the cut-off momentum t m i n in the scattering subprocesses. 
Therefore, the contribution of the charm excitation mechanism in the hadron collisions is, as 
minimum, order of magnitude less than that obtained in [8] and cannot explain large yield of 
charmed particles, observed in ISR experiments. Thus, the contradiction between ISR data 
and data of fixed target experiments is not avoided. 

As it follows from the results of present work, it is not necessary to introduce any visible 
admixture of the non-perturbative charm quarks in the proton to describe EMC data, even when 
using "soft" gluon distribution. The excess of the EMC data over the photon-gluon predictions 
at large Q 2 can be completely explained by the contribution of the hadronic component of the 
photon. In any case, the considered mechanism decreases significantly the discrepancy between 
EMC data and PGF predictions and, respectively, the admixture of Fock states with charmed 
quarks in the proton. 

On the other hand, we have introduced a non-perturbative charmed component in the pho- 
ton wave function. The reason, why the photon can have visible non-perturbative component 
of c quarks in difference to the proton can be understood from result of paper [24] , which states 
that the perturbative regime in the photon starts at Q 2 pa 5.5 GeV 2 , and for the proton the 
DGLAP equations give good description of the evolution of the structure function starting al- 
ready from Q 2 pa (1 — 2) GeV 2 . Thus, there is no necessity to introduce the non-perturbative 
charmed states in the proton and the charm contribution to the proton structure function can 
be described by pQCD. 
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4 Appendix 

Expressions for total and differential cross sections of c quark scattering at lowest order of a s 
[30] (we write out these formulae as they are given in [8]). 



da(qq^> cc) 

~dl 



An „,2{a 2\ (m 2 -i) 2 + {m 2 .-u) 2 +2m 2 s) 



a(qq -> cc) = |^a 2 (4m 2 ) (s + 2mc 2 ) 
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s th = 4 m; 



dcr(gg^cc) 
di 



= i§2a s 2 (4m;D 
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i« s 2 (4m c 2 ) 



7 + 
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16(1— to 2 )- 



tmin) 



32(s-mg)(mg-M) 
£ 2 



+ 



64 (g~TO 2 J(m 2 -n)+2TO 2 (s+m;j) 64 {s-ml){m^-u)+2ml(u+m 2 c ) 
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i(s-m 2 ) 



- 16 



(J— to 2 ) (to 2 — u)— to 2 (s— u) 
t(m 2 — ?1) 



(s-to 2 ) 2 «s( *min) 



^ + 9 (s-to 2 ) J (*min 



tmax) 



9 f 2 ■ — f 2 
I vain max 

' 9 s-to 2 



+ 2(s + m2)log|^ 



. 4 s 2 -6to 2 s+6to^ i s-TO 2 +£ m m 



s— m; 



s-to 2 +£„ 



+ 2(«- m ») (r 5 --^) 

\ L max L min / 



l 



9 c \ s— m 2 +<„ 



'max 



= - max (s - mj? + f min , (s - rrfy/s) , 



(to enforce u — < £ mm ) ; 



2 IX / 2 ^ 1 X2 \ • r X 
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171X71 J 



if 



mm 



16 



(a) (b) 




fl 



(c) (d) 



Figure 1: Processes considered in present paper and giving a contribution in structure function 
of the charm leptoproduction: (a) — the lepton scattering off c quark (the radiative corrections 
are not depicted); (b) — photon-gluon fusion; (c), (d), (e) — the scattering of c quark from 
the photon on gluon and quark from the proton, respectively. 
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Figure 2: A comparison of the result obtained in present work with experimental data of EMC. 
Solid line - sum of all contributions, dashed line - the contribution of the hadronic component 
of the photon, dotted line - PGF, dot-dashed line - the scattering of the point-like photon on 
c quark from the proton ("intrinsic charm"). For the gluon distribution the parametrization 
MRS(G) was used. The results are given for the following values of the parameters: i min = 
-3.01 GeV 2 , N IC = 0.20% and mass of c quark m c = 1.51 GeV. 
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